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When  a  lightly-damped  osclUator  Is  excited  by  s  stationary  random 
process  the  response  has  the  aipearance  of  an  oscillation  vlth  a  slow, 
random  fluctuation  In  aiqplltude.  The  envelope  of  the  oscillation  has 
been  st\idled  by  Rice  |l^  for  the  case  of  stationary  Oausslan  response. 
When  Gaussian  excitation  Is  applied  to  a  nonlinear  oscillator  the 
response  Is  generally  non-Gausslan .  In  the  non-Gaussian  case  the  Rice 
technique  for  defining  an  envelope  Is  no  longer  directly  applicable. 

It  is  not  clear  vhether  there  la  any  possibility  of  extending  the  Rice 
approach  to  obtain  mean:,'jjgfUl  results  In  the  non-Gausslan  case.  The 
present  author  has  given  two  alternative  techniques  [2J  which  permit 
evaluation  of  the  first-order  probability  distribution  for  the  response 
envelope  of  a  nonlinear  slngle-degree-of- freedom  oscillator  with 
Gaussian  excitation.  The  purpose  of  this  note  Is  to  show  that  for  a 
general  class  of  nonlinear  restoring  forces  these  two  first-order 
distributions  are  Identical.  Fuztbenoore,  for  linear  systems  this 


distribution  Is  Identical  with  that  given  by  the  Rice  approach.  Our 


ABSTRACT 


The  Envelope  of  Random  Vibration  of  a 
Lightly-Damped  Honlinear  Oacillator 


by 

Stephen  H.  Crandall 
Professor  of  Mechanical  Engineering 
Massachusetts  Institute  of  Technology 

There  are  several  ways  in  which  the  envelope  of  the  stationary 
random  response  of  a  lightly-dan^d  oscillator  can  he  characterized. 
Two  characterizations  which  lead  to  the  same  first-order  prohahility 
distributions  for  systems  with  conservative  nonlinear  restoring 
forces  are  described.  In  such  systems  the  envelope  distributions 
are  in  general  different  from  the  distributions  of  peak  amplitudes. 
The  Duffing  system  is  used  to  illustrate  these  general  results. 

For  the  special  case  of  a  Gaussian  oscillator  an  envelope 
characterization  has  been  given  by  Rice,  It  is  shown  that  In  this 
case  the  above  envelope  characterization  has  the  same  first-order 
distribution  as  that  given  by  Rice,  but  that  the  second-order 
probability  distributions  are,  in  general,  different. 
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enrelope  definition  is  never- the- lees  not  completely  identical  with 
that  of  Rice  even  for  narrow-hand  Gaussian  responses  because  the  second- 
order  probability  distributions  for  an  id'-al  band-pass  filter  differ 
by  tenas  proportional  to  the  square  of  the  bandwidth. 


1.  Envelope  Obtained  from  Peedi  Statistics 

We  consider  nonlinear  oscillators  described  by  the  equation 


;V  I- fi  x  ■+  ^(x)  =.  (ic) 


(1) 


where  ^  is  a  constant  linear  damping  coefficient,  g(x)  is  the  nonlinear 
restoring  force  function  emd  f(t)  is  the  excitation  vdiich  ve  will  take 
to  be  a  stationary  Gaussian  process  with  uniform  spectral  density 
(acceleration  sc^uared  per  unit  of  circular  frequency)  over  all  frequencies 
from  4^  =  -OOto  We  will  assume  that  g(x)  is  an  odd  function 

such  that  xg(x)  is  positive  definite.  The  integral 


(2) 


is  thus  a  positive  definite  even  flinction.  It  is  further  appropriate 
to  apsujrsc  that  g(x)  remains  sufficiently  large  for  large  x  that 


ri  a/ 


<  oo 


(3) 


where  ^  is  a  real  constant.  The  condition  (3)  insures  that  the 
stationary  response  x(t)  of  (l)  has  a  finite  mean  sqxiare. 

The  Joint  probability  distribution  density  yi7(x,v)  of  x  and  v  =  x  for 
the  stationary  response  of  (1)  can  be  obtained  as  a  solution  of  the 


-3- 


Fokker-Planck  equation  |_3j  •  We  find 


h  I 


K'V)  -  c 


(»*) 


Where  C  i«  the  normalization  constant  and 

tS^ 


I' 


(5) 


is  the  mean  sq^uare  value  of  v(t)  in  the  stationary  response.  The  expected 
number  of  crossings  of  the  level  x  =  a  with  positive  slope  per  unit 
time  is  [l] 

oo 

7^^  ~  ^  jo  iCl,rj)  (6) 


Inserting  (4)  into  (6)  yields 


hi/ 


e 


G(^) 


-z 


(7) 


where  is  the  exi)ected  number  of  zero-crossings  with  positive  slope 
per  unit  time.  The  statistical  parameter  (7)  can  also  be  considered  to 
represent  a  probability  distribution  by  the  following  heuri«ti<'  argument. 
In  a  narrow-band  process  the  response  consists  essentially  of  recognizable 
cycles  with  one  positive  peak  and  one  negative  peak  peycycle.  The 
parameter  counts  the  average  number  of  cycles  per  unit  time  while 

the  parameter  counts  the  ave:rage  number  of  cycles  per  xinlt  time  with 
peaks  above  the  level  x  =  a.  Thus,  on  the  average,  the  fraction  of  cycles 


with  peaks  higher  than  x  s  a  Is 


p  (  jo&ok  >  d) 


The  probability  density  of  peaks  is 


dbCu 


e 


-tr 


substituting  from  (7)  and  (2).  It  is  to  be  emphasized  that  (8)  and  (9) 
represent  the  distribution  of  peaks  within  the  population  of  all 
positive  peaks. 

An  envelope  distribution  can  be  obtained  from  (9)  by  assuming 
that  the  envelope  is  a  smooth  gradual  curve  Joining  the  peakt.  :..ad  that 
on  the  average  the  time  the  envelope  spends  between  a  and  a  +  ^  is  Just 
the  number  of  such  peedts  multiplied  by  TT  (a)  which  is  the  undamped 
period  of  a  free  vibration  of  amplitude  a.  The  fraction  of  time  so 
spent  is 

(b  (  a)  ia  =  /)  W)  tI*  djx 

6  ^ 


The  undamped  period  (a)  can  be  obtained  by  Integrating  (1)  with  0 

and  f(t)  =0. 

a 

\/T[G(Ci)  - 
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Inserting  (9)  and  (11)  Into  (10)  we  have  the  first-order  probability 
density  for  an  envelope  of  the  response  process  defined  by  (l). 


Mm 


(V 


CL 


d'x 


\/  -  Gl'xj] 


(12) 


Evaluation  of  (12)  for  a  particular  nonlinear  fliactlon  g(x)  may  be 
difficult.  In  addition  to  the  operations  explicitly  Indicated  in  (12) 
the  parameter  d  req^uirec  evaluation  of  the  normalization  constant  in 
(V)  as  follows 


(13) 


The  arguments  leading  to  (12)  are  not  rigorous.  Fuithermore,  the 
envelope  so  defined  is  insufficiently  detailed  to  permit  extension 
beyond  the  first-order  probability  distribution  eind  its  accompanying 
first-order  statistics.  We  next  consider  a  more  definitely  defined 
envelope  \rtiich  turns  out  to  have  the  san^  first-order  distribution 
as  (12)  . 


3*  Envelope  defined  by  total  energy 

Along  with  the  random  processes  x(t)  and  v(t)  determined  by  (l) 
let  us  consider  the  envelope  process  a(t)  vhere 


G(a)  ---  Qi^) 


2 


(14) 


-6- 


The  interpretation  here  is  tL\at  a(t)  ia  the  amplitude  which  would  be 
attained  if  the  sum  of  the  instcmtaneous  kinetic  and  potential  energies 
were  converted  entirely  into  potential  energy;  i.e.,  a(t)  is  the  peak 
amplitude  \dilch  would  be  obtained  in  a  hypothetical  free  undamped 
oBcillatlon  of  the  system  if  the  initial  conditions  for  the  free 
oscillation  were  the  cxxrrent  values  of  x(t)  and  v(t) .  The  first-order 
probability  distribution  for  a(t)  then  follows  from  the  Joint  distribution 
of  X  and  v.  Taking  advantage  of  symmetry  we  have 


P  [emveio/Dt  <  (Z) 


4 


dx 


/  d'lr 


(15) 


as  the  ensemble  fraction  of  cases  for  which  the  combination  of  x  and  v 
makes  the  envelope  defined  by  (lU)  smaller  than  a  given  value  of  a. 

The  probability  density  for  the  envelope  is 


CU 


3(0.}  dx 


(16) 


Inserting  the  Joint  density  (U)  leads  to 


which  is  identical  with  (12)  if  we  take  into  account  the  relation  (13) 
between  the  normalization  constant  C  and  the  average  frequency  •  Thus 
for  arbitrary  nonlinear  restoring  forces  of  the  type  indicated  the  heuristic 
envelope  of  Section  1  has  the  same  first-order  probability  distribution 
as  the  envelope  based  on  the  energy  relation  (1^). 
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3.  ExaBg>lg  -  Duffing  System 

We  consider  the  particular  case  of  (l)  In  which 


<3(x)  =  ( i  ^  f  } 


where  It  Is  possible  to  evaluate  the  general  expressions  obtained  above. 
A  convenient  normalisation  Is  provided  by  the  parameter  ( where 

tSc 


C  - 


Is  the  mean  square  value  of  In  the  limiting  linear  case  \diere  L  -9/0 
Setting 


I 


f  -  S 

^  L 


the  governing  eq,\iatlon  (1)  takes  the  form 

J  (ij  ^  ^  = 


The  free  undamped  osclUationB  of  angjlitude  yjn^x  *  ®  ^  described 

In  terms  of  the  Jacobleui  elliptic  cosine  Sanction 


^  -  a  cn  ^ 


where  m  is  the  moduhis  of  the  elliptic  function 


m  - 


and  the  parameter  (J  Is  given  by 


6  'CC 


2L1  ^  e^a^j 


2  2-  /  £^2 


(18) 


(19) 


(20) 


(21) 


(22) 


(23) 


(24) 
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The  period  'X^{a)  is 


where  K(m)  is  the  complete  elliptic  Integral  which  giviss  the  real  q^uarter 
period  of  the  elliptic  ftinction  [kj  . 

The  joint  probability  b(y,y)  for  this  case  follows  from  (4) 


Ky^y)  = 


{St) 


where  Integral 


provides  the  correct  no rmsli nation.  The  probability  density  for  the  peaks 
(9)  is  obtained  without  evaluating  (Z/) 


a 


a(  / 


^  a 


e 


m 


and  is  plotted  in  Fig.  1  for  several  values  of  the  nonlinearity  parameter  £,  . 
This  distribution  was  recently  obtained  by  Lyon  [5J  using  essentially  the 
same  procedure  as  ttet  described  in  Section  1.  In  it  was  additionally 
suggested  that  (Z6)  might  also  be  interpreted  as  the  density  function  of  an 
envelope.  Since  this  suggestion  takes  no  account  of  the  tendency  of  ^cycles” 


with  different  peak  anplitiides  to  have  different  periods  we  doubt  that  it 


portrays  an  envelope  as  effectively  as  the  distribution 


^  X(.a)  t>^((i) 

Xl^)\I  i+€^aJ^ 

which  follows  from  (10)  by  using  (2’5)  and  the  relation 

_  'A 

K^) 

which  is  obtained  from  (13)  •  The  distribution  (29)  can  also  be  obtained 
directly  from  (17)  when  the  Integral  is  recognized  as  the  complete 
elliptic  integral  of  (2*5),  Fig.  2  shows  plots  of  (2LP)  for  different 
values  of  C  . 

For  the  Duifing  oscillator  there  is  not  a  great  difference  between 
{2S)  and  (29)  until  £  is  quite  large-  In  Fig.  3  there  is  a  comparison 
between  the  distributions  of  (£5)  aind  (29)  for  £  =  2*  When  £3  0  the 
system  {2f)  is  linear,  the  undamped  free  vibration  period  is  Independent 
of  anqplitude  and  the  distributions  (2S)  and  (25)  both  reduce  to  the 
Rayellgh  distribution 

lD(a)  ~  d  ^ 

which  is  the  known  amplitude  distribution  for  narrow-band  Gaussian 
processes  and  serves  equally  "well  to  ucbcrlbe  T.he  peaks  or  zne  envelope. 

Thfi  first-order  distribution  of  the  response  ^  is  of  some  auxiliary 
interest  for  this  system.  It  follows  from  (2^)  by  integration  over  all 


values  of  y. 
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e 


This  is  shown  in  Figure  h  for  several  values  of  the  nonlinearity  parameter 

€  • 


ii-.  Second-order  Envelope  Distribution 

The  energy-letsed  enveloi>e  definition  (l4)  gives  promise  of  providing 
a  means  for  obtaining  second-order  probabilities  or  second-order  statistics 
for  the  envelope  of  the  response  of  (l)  to  Gaussieui  excitation.  The  joint 
probability  that  at  time  jt  the  envelope  will  oe  less  than  a^  and  that  at 
time  t  +~l^  the  envelope  will  be  less  than  is  given  by 


P  =  \ 


dv-^ 


r 


J 


J 


\/^CQo.^yG(^z)j 


(>23) 


as  a  direct  extension  of  (15)-  The  density  p(a^,a2)  then  foUovs  from 
differentiation  of  (33)  •  Our  work  along  this  line  has  so  far  been 
limited  to  cases  where  p(x^,Vq^,X2,V2)  has  been  Gaussian.  In  the  Gaussian 
case  we  will  show  that  the  second-order  density  p(aijtap)  so  obtained  is 
not  identical  with  the  corresponding  density  which  foDJLovs  from  Rice  *8 
envelope  definition. 

When  the  system  (1)  is  linear  the  response  to  Gaussian  excitation 
is  also  Gaussian.  In  this  case  the  enveloi)e  definition  (l4)  ceui  be  altered 
to  make  it  depend  solely  on  the  response  process;  1.©.^  independent  of  the 
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oscillator  parameters.  Thus  for  a  narrow-lxand  Gaussiem  processes  x(t) 
the  energy-based  envelope  a(t)  is  given  by 


CL  - 


% 


■i- 


{%) 


2. 


(34.) 


)  + 

where  =  2  is  the  expected  circular  frequency  of  the  process 

and  V  ■  X.  The  result  (34)  follows  from  setting  G(x)  x'^  in  (l4) 

The  first-order  distribution  of  a(t)  is  that  of  Rayleigh  as  noted  in  the 
previous  section.  The  second  order  density  follows  from  (35)  but  because 
of  (34)  a  great  sin^Jllflcatlon  is  obtained  b>r  transforming  to  polar 
coordinates .  We  find 


r 


?Tr 


ijfr 


In  order  to  evaluate  this  it  is  necessary  to  know  the  Joint  density 
p(xq^,v^,X2,V2)  ^ich  in  turn  depends  on  the  autocorrelation  function 


R(2-;  =  EiiXjxj 


and  its  first  two  derivatives.  We  shall  limit  our  further  discussion  to 
a  particular  choice  for  {3i>) . 


5“  Example -Ideal  Band  Pass  Filter 

Let  the  spectral  density  of  x(t)  be  zero  except  in  bands  of  width  2<7C 
centered  on  +  6l^ where  the  density  is  constant.  The  correspondinK  auto¬ 
correlation  function  is 


R  [T)  =  cosd„z 


(37) 
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where  (f  ^  =  E  Is  the  mean  square  of  the  process.  We  introduce  the 


follovlng  notation 


S/1  9^ 


=  3 


(3») 


/  “  10-  CJ, 


The  functions  fj^  and  fg  are  shown  in  Figure  5-  The  parameter  ^  is  a 
measure  of  the  process  bandwidth.  In  terms  of  these  we  find  successively 

'R(Z)  =  (?■'■  l;(o<Z)C0JC4)^Z 

■  P  ix)  -  (T  -h  X  fiUT)^r>fX  ^ 

{  ^ 


h 


a\/in 


icfo  trr>^rLS 


The  joint  density  of  four  Gaussian  variables  y^,  yg,  y^  and 


\1 


(4^) 


vrtiere  iB  th<5  deteminant  and  the  Aj  are  the  cofactors  of  the 
covariance  matrix  for  the  y^.  If  we  identify  the  with  x^,  v^,  x^ 
respectively  the  covariance  matrix  is 


“7 


R(0) 

0 

R(r) 

B-ra 

0 

-R"(0) 

-R’(r) 

-R"m 

B(T) 

-B’(r) 

-R(0) 

0 

B'(T) 

-B"(r) 

0 

-R"(0) 

and  V 

2 
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In  principle  all  that  remains  is  to  insert  (3®)  into  {4l  evaliaate  ^o) 
and  then  obtain  the  second-order  probability  distribution  density  using 
{35)  •  have,  however,  been  unable  to  complete  this  program  without 
approximation;  but  if  we  neglect  terms  which  are  0(  ^  )  in  (55)  the 
calculation  is  rather  simple.  For  the  distribution  {^)  we  find 


6r/ 


-/V 


2(r*4,'-(r-  i'j 


{4a) 


with 


/[/  —  ^  ~  2  (. 

\Tt^ 

-f*  2  ^  5  m  C  "<'7] 


Here  the  argument  of  f-^  is  as  in  (33)- 

A  final  preliminary  before  evaluating  (33)  is  the  determination  of 
the  expected  frec^uency 


- 

R(o) 


( 


\44i 


P  \ 

Here  again  we  have  neglected  the  v  term  and  taken  to  be  identical 

with  .  The  Integration  of  (32)  then  becomes  straightforward  and 

we  obtain  for  the  joint  density 


hM  a)  -  — T 


cUfa^-C 


I 

o  L  r  o-i 


(43J 


which  iB  Identical  with  that  obtained  by  Rice  jlj  with  a  dj^fe  •  nt  en/elope 
definition.  Since  we  have  had  to  make  0( ^  )  approximations  whereas  no 
such  approximation  was  required  in  jlj  there  is  a  Buggestion  tha^  the 
two  envelopes  are  not  the  same . 

In  order  to  establish  this  point  we  turn  to  a  simple  second--  rder 
envelope  statistic  which  can  be  evaluated  exactly  for  both  envelope 

definitions.  Using  Rice's  definition  and  hence  (.^5) 

ao  pcC> 

=  (  cC  dA,  ( 

<Jo 

^  4-  ^  1 

where  the  integration  is  performed  along  the  lines  indicated  in  f^sj  and 
use  is  made  of  jjJ,  j^J  and  jyj. 

Using  the  definition  (34)  have 


i«]  =  eCsYJ 


and  since  x-i ,  Xp,  vi  and  Vp  are  Jointly  Gaussian 


EC-KfCj  - 

£  - 


(T*^  +  2 

cj*<r'*  +  2  [  R'dr)]- 


SO  thac  using  (59)  Permits  us  to  evaluate  (47)  exactly.  The  result  is 
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fairly  complicated  but  It  can  be  put  la  the  following  form 

EroiVJ  =  4-r^  -t 


with 


fyz)  -  jcscos^cz^  ~ 


which  indicates  that  there  definitely  is  an  0{y)  discrepancy  between 

r  2  2  n 

the  values  of  E  la a2  ^  obtained  from  Rice’s  envelope  definition  and 
from  the  energy-based  envelope  (^34)  •  Finally  it  should  be  noted  that 
this  discrepancy  applies  only  to  the  special  ideal  band-pass  process  of 
(37)  -  For  other  narrow-band  processes  the  discrepancy  may  not  necessarily 
be  pi"oportional  to  the  scLuare  of  the  bandwidth. 
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CAPTIOHS  FOR  FIGURES 


1.  Probalsillty  density  distribution  of  peaks  for  Duffing  oscillator 
subjected  to  uhite-noise  excitation  for  various  values  of  the 
nonlinearity  parameter. 

2.  Probability  density  distribution  of  an  envelope  for  the  Duffing 
oscillator. 

3*  Comparison  of  probability  density  distributions  of  the  peaks  {t8) 
and  of  the  envelope  (23)  for  the  Duffing  oscillator  with  €  »  2. 

U.  Probability  density  distribution  for  Duffing  systeiz  response  process. 

5.  The  functions  of  and  ^2^^^  defined  in  (33)  and  used  to 

describe  the  correlations  for  the  ideal  band-pass  filter  response. 


FIGURE  3 
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